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ABSTRACT 


% 


SRL  has  designed,  fabricated,  and  successfully  tested  a  pulsed  C02  laser  system.  <  Laser 
pulses  have  been  generated  with  measured  output  energy  of  600  mJ,  pulse  length  40  /usee,  and 
conversion  efficiency  of  14%  in  a  3:2:1  (He:N2:C02)  gas  mixture.  The  laser  operates  at  a  total 
gas  pressure  between  5  and  40  Torr,  and  utilizes  the  COLD-I  all-solid-state  C02  laser  driver 
developed  by  SRL.  Typical  current,  voltage,  and  laser  output  traces  are  presented  and  compared 
to  numerically  simulated  curves,  showing  reasonable  agreement. 

The  theoretical  model  of  C02  laser  discharges  developed  by  SRL'has  been  extended  and 
utilized  to  analyze  the  steady  state  solutions  and  stability  criteria  for  three  regimes  of  operation: 
the  stiff  voltage  case,  the  stiff  current  case,  and  the  finite  external  impedance  case.  Numerical 
solutions  are  presented  for  the  temporal  evolution  of  the  electron  (plasma)  density,  the  metastable 
density,  and  the  electric  field  under  various  conditions.  The  theory  of  RF  discharge  stabilization 
has  been  extended  to  two  spatial  dimensions,  and  numerical  results  are  presented  for  the  stiff 
current  and  the  stiff  voltage  cases. 

In  a  repetitively  pulse  C02  laser,  the  lasing  gas  must  be  cooled  between  pulses  to  preserve 
efficiency.  Conduction  cooling  has  been  examined  by  SRL  as  an  advanced  alternative  to  the 
traditional  flow  loop  approach.  This  concept  minimizes  the  weight  of  the  laser  head.  A  thermal 
analysis  of  a  pulsed  C02  laser  head  with  conduction  cooling  is  presented,  along  with  a  discussion 

f 

of  beam  quality  effects  due  to  acoustic  and  thermal  processes. 
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SECTION  1 


1.0  RESULTS  FROM  THE  SRL  C02  LASER 

During  the  period  of  this  report,  SRL  has  completed  the  design,  fabrication,  and  initial 
testing  of  a  1  meter  conduction  cooled  C02  laser  head.  In  this  section,  the  laser  head  and  the 
all-solid-state  C02  laser  driver  are  described,  and  performance  results  are  compared  to  numerical 
simulations. 

1.1  Laser  Design 

An  assembly  drawing  of  the  C02  laser  head  is  shown  in  Fig.  1.1.  The  discharge  volume 
is  1  meter  long,  and  the  anode-cathode  gap  is  13.5  cm.  The  discharge  width  is  limited  to  2.0 
cm  by  glass  side  walls  within  the  laser  head.  This  spacing  is  designed  to  enable  the  removal  of 
heat  from  the  laser  gas  under  repetitively  pulsed  operation  by  conduction  cooling,  as  discussed 
in  Section  3.  Our  calculations  predict  that  operating  the  laser  at  a  10  pps  repetition  rate  will 
result  in  a  steady  state  temperature  rise  of  75°C,  which  will  not  degrade  the  laser  operation. 

The  discharge  electrodes  have  smoothly  tapered  edges  along  the  sides  and  at  the  ends  to 
prevent  discharge  instabilities  and  arcing.  Embedded  within  the  cathode  (but  electrically  isolated) 
are  numerous  preionization  sites  consisting  of  short  wire  pins  to  which  a  prepulse  is  applied  in 
order  to  uniformly  initiate  the  discharge.  To  date,  the  laser  has  been  operated  with  22  such 
pre ionization  sites.  Typical  discharges  in  10-40  Torr  of  the  He  or  He/N2/C02  gas  mixtures  show 
an  expanding  discharge  emanating  uniformly  from  the  immediate  region  of  the  wire  pins  towards 
the  anode. 

The  electrode  profile  is  a  modified  contour  designed  to  eliminate  arcing  problems  originally 
seen  near  the  electrode  edges.  The  modified  profile  yields  a  stable  discharge.  Both  contours  are 
shown  in  Fig.  1.2,  along  with  a  schematic  of  the  preionization  pins.  Each  preionization  pin  is 
individually  ballasted  and  then  connected  to  one  common  capacitor.  It  has  been  demonstrated 
that  all  the  pins  fire  simultaneously,  as  long  as  the  ballast  resistors  exceed  20  Ohms. 

In  addition  to  the  laser  itself,  the  system  includes  the  pulsed  power  supply,  a  gas  handling 
system  including  a  mixing  tank,  and  a  CAMAC  data  acquisition  system.  These  systems  were 
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completed  and  tested  concurrently  with  the  laser.  The  power  system  and  the  data  system  are 
described  below. 

1.2  Pulsed  C02  Laser  Driver 

The  design  of  the  COLD-I  (C02  Laser  Driver)  pulsed  power  supply  has  been  presented  in 
detail  in  a  previous  report.  The  COLD-I  pulser  is  a  compact,  lightweight  all-solid-state  device 
employing  SCR  commutators  combined  with  nonlinear  magnetic  pulse  compression.  During  this 
reporting  period,  the  pulser  has  been  fabricated,  tested,  and  incorporated  into  the  laser  discharge 
circuit.  A  photograph  of  the  COLD-I  pulser  with  a  small  laser  head  is  shown  in  Fig.  1.3.  The 
pulser  operated  properly,  as  expected.  The  I/V  characteristics  of  the  main  discharge  pulse  are 
discussed  below.  A  separate  preionizing  pulse  is  applied  to  the  22  wire  pins  immediately  prior 
to  the  main  pulse. 

1.3  Diagnostics  and  Data  System 

Three  diagnostic  systems  were  installed  to  measure  the  discharge  current,  the  discharge  volt¬ 
age,  and  the  laser  optical  output  power  as  functions  of  time.  The  current  pulse  is  monitored 
by  a  Pierson  wide  band  current  transformer.  The  anode-cathode  voltage  is  detected  by  a  fiber 
optic  voltage  probe  which  is  immune  to  troublesome  inductive  voltage  swings  on  the  nominally 
grounded  anode.  This  probe,  shown  in  Fig.  1.4,  consists  of  two  LED  fiber  optic  transmit¬ 
ter/receiver  loops.  The  resistor  chain  limits  the  current  through  the  LED’s  to  approximately  40 
mA  at  the  maximum  voltage  on  the  laser  head.  One  LED  will  send  out  its  light  signal  when 
the  voltage  across  the  laser  head  is  positive  and  the  other  will  send  out  its  signal  when  the 
voltage  difference  is  negative.  There  is  a  very  small  dead  zone  near  zero  voltage  difference 
since  the  diodes  require  a  volt  or  so  of  forward  bias  before  any  conduction  occurs,  but  this  is  a 
tiny  fraction  of  the  kilovolt  swings  being  measured.  The  diagnostic  has  an  unambiguous  ground 
reference  and  noise  is  reduced  since  there  is  no  direct  electrical  connection  between  the  LED 
voltage  probes  and  the  data  acquisition  system.  This  high  voltage  probe  was  shown  to  be  linear 
up  to  at  least  15  kV  during  calibration  studies  with  a  pulse  generator.  Reliable  current  and 
voltage  measurements  are  obtained  since  both  the  voltage  and  current  diagnostics  are  decoupled 
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from  the  floating  ground  on  the  laser  pulsed  power  supply. 

Finally,  the  optical  laser  output  pulseshape  and  energy  are  measured  by  a  photon-drag  de¬ 
tector  with  1  nsec  response  time  at  the  C02  wavelength  and  a  pyroelectric  energy  meter,  re¬ 
spectively.  All  diagnostic  outputs  are  recorded  by  a  data  acquisition  system  assembled  from 
modular  CAM  AC  digitizers.  Three  100  megasample  sec-1  digitizers  have  been  installed  to  read 
voltage,  current  and  the  laser  pulseshape.  Three  32  megasample  sec-1  digitizers  are  available 
for  monitoring  the  output  of  the  pyroelectric  energy  meter  and  other  test  point  voltages  in  the 
pulsed  power  system.  A  data  acquisition  program  was  written  to  program  the  digitizers,  acquire 
the  voltage  traces  and  store  the  data  in  standardized  file  structures.  The  program  provides  for 
automatically  generated,  unique  file  names  and  appends  descriptions  of  each  test  run  to  the  indi¬ 
vidual  data  files.  Provision  has  been  made  for  recording  2-D  images  of  beam  profiles  or  of  arc 
structure  with  a  CCD  camera  if  necessary. 

1.4  Discharge  Characteristics 

After  initial  testing  of  the  COLD-I  pulsed  power  unit  with  CuS04  resistive  loads  ranging 
from  4  to  50  Ohms,  the  pulser  was  used  to  drive  gas  discharges  in  the  laser  head.  With  a 
3:2:1  (He:N2:C02)  gas  mixture,  the  response  of  the  gaseous  discharge  to  the  pulser  has  been 
investigated  intensively.  Total  gas  pressures  from  5  to  40  Torr  have  been  studied  The  input 
voltage  to  the  pulser  was  varied  from  500  to  1100  V;  the  discharge  characteristics  were  found 
to  be  optimized  at  the  higher  input  voltages. 

Voltage  and  current  traces  from  a  20  Torr  discharge  at  1000  V  input  are  shown  in  Fig.  1.5. 
The  peak  output  power  to  the  discharge  is  12  MW  with  an  impedance  of  4  Ohms.  Voltage 
ringing  during  the  first  //sec  is  followed  by  a  self-sustained  discharge  lasting  an  additional  2-3 
//sec. 

A  35  Torr  discharge  is  shown  in  Fig.  1.6,  also  taken  with  1000  V  input  to  the  COLD-I 
pulser.  A  numerical  simulation  of  this  discharge  was  performed  to  model  the  electron  density 
and  total  current  as  a  function  of  time.  The  simulation  is  based  on  two  rate  equations: 
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CO  2  Laser  Discharge  Characteristics 
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Figure  1.6  Voltage  and  current  traces  from  35  Torr  discharge, 
including  current  predicted  by  numerical  simulation. 
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where  ne  is  the  electron  density,  na  the  density  of  N2  in  the  ground  state,  nc  the  density  of 
CO2  in  the  ground  state,  nm  the  density  of  C02  molecules  in  an  electronic  metastable  state,  v 
the  ionization  rate  of  the  metastable  C02,  ua  the  ionization  rate  form  the  ground  state  of  N2,  a 
the  ion-electron  recombination  rate,  (3  the  electron  attachment  rate,  7  the  collisional  excitation 
rate  for  the  metastable  C02,  and  r  is  the  lifetime  of  the  metastables.  From  the  electron  density 
and  voltage  across  the  electrode  one  obtains  the  current  through  the  electrode.  The  results  of 
this  simulation  are  also  shown  in  Fig.  1.6.  The  close  agreement  between  the  simulation  and 
experimental  results  suggest  that  Eqs.  1.1  and  1.2  describe  the  discharge  process  very  well. 

The  discharge  was  stable  and  no  arcs  were  observed  along  the  insulators.  An  open  shutter 
photograph  of  the  20  Torr  discharge  is  shown  in  Fig.  1.7.  From  the  photograph  it  is  clear  that  the 
discharge  current  density  is  larger  in  the  vicinity  of  the  preionizers.  The  redesigned  electrodes, 
which  increase  the  number  of  preionization  pins  from  22  to  120,  will  improve  the  uniformity  of 
the  discharge. 

1.5  Laser  Performance  Results 

The  C02  laser  cavity  consists  of  a  molybdenum  back  mirror  (uncoated)  having  a  10  meter 
radius  of  curvature  and  a  ZnS  50%  transmission  output  coupler  (AR  coated).  Typical  laser 
outputs  are  shown  and  compared  to  numerical  simulations  for  the  20  Torr  and  35  Torr  discharges 
discussed  above  in  Figs.  1.8  and  1.9,  respectively. 

At  20  Torr  (Fig.  1.8),  the  observed  laser  output  decays  exponentially  with  a  time  constant  of 
28  fi sec,  which  corresponds  to  the  collisional  energy  transfer  time  between  N2  and  C02.  This  is 
in  good  agreement  with  the  numerical  simulation  for  that  pressure,  which  exhibits  a  decay  time 
of  25  n  sec. 

At  35  Torr  (Fig.  1.9),  the  laser  output  builds  up  on  a  2  /zsec  time  scale  after  the  initial  0.5 
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CONDUCTION  COOLED  C02  LASER  DISCHARGE 


Active  Volume 

Voltage 

Electrode  Spacing 
Current 

Discharge  Pulse  Length 
Laser  Pulse  Length 
Overall  Efficiency 
Total  Gas  Pressure 


2.0  x  13.5  x  100  cm'' 
5  kY 


13.5  cm 
S00  Amps 
3  fis 
40  //s 
14'/ 


30  Ton 


Figure  1.7 
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3:2:1(He:N2:C02  )  35Torr 
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fiscc  spike,  and  has  an  exponential  decay  time  of  21  fiscc  due  to  the  faster  collision  rate.  The 
tail  of  the  pulse  is  over  40  //sec  long.  This  agrees  very  well  with  the  simulated  pulse  shape, 
which  has  a  decay  time  of  16.4  /xsec.  The  simulation  exhibits  more  ringing  of  the  pulse  after 
the  initial  spike  than  is  detected  experimentally.  The  maximum  total  laser  output  energy  is  600 
mJ,  corresponding  to  a  conversion  efficiency  of  14%. 
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SECTION  2 


2.0  LASER  DISCHARGE  STABILIZATION 
2.1  Introduction 

In  previous  reports,  SRL  has  presented  a  stability  model  for  one-dimensional  DC  volumetric 
discharges.  Stability  was  found  to  depend  critically  on  the  electron  attachment  rate  constant 
(/?)  for  several  operating  regimes.  In  this  section,  those  results  are  summarized  and  extended 
in  terms  of  the  critical  electric  field  required  for  stable  operation  in  each  regime.  Kinetic 
rate  coefficients  previously  calculated  and  reported  for  a  one  atmosphere  total  pressure,  3:2:1 
(He:N2:C02)  mixture  are  used  in  numerical  calculations  presented  below. 

In  addition  to  the  DC  stability  model,  SRL  has  developed  a  two-dimensional  model  for  RF 
discharges  which  demonstrates  the  possibility  of  inductive  stabilization  of  the  discharge  against 
streamer  formation.  This  model,  reported  previously,  has  now  been  incorporated  into  a  numerical 
code  and  used  to  predict  discharge  behavior  under  conditions  of  RF  electric  fields.  These  new 
results  are  presented  in  Section  2.5  below. 

The  stability  of  high  pressure  discharges  has  been  a  subject  of  major  concern  in  the  high 
power  laser  community  over  the  last  fifteen  years.  Most  of  the  discharge  stability  research, 
especially  for  C02  laser  discharges,  concentrated  on  thermal  instabilities  and  their  effect  on  the 
constriction  and  arcing  of  the  discharge.  Vibrational  stabilization  methods  have  been  proposed 
to  suppress  these  instabilities. 

The  thermal  instabilities,  by  their  nature,  can  only  develop  on  acoustic  (millisecond)  time 
scales.  However,  in  many  cases  arcing  occurs  on  a  much  faster  (micro-second  and  submicrosec¬ 
ond)  scale.  The  phenomenon  is  characteristic  of  systems  where  the  electron  loss  by  attachment 
is  present  and  the  two-step  electron  impact  ionization  is  dominant.  Such  a  model  success¬ 
fully  explained  the  experimentally  observed  temporally  rapid  instability  scales  in  the  KrF  laser 
discharges. 

In  addition  to  the  KrF  laser  discharge  case,  subacoustic  arcing  time  scales  have  been  also 
observed  in  the  high  power  C02  laser  experiments.  An  example  of  such  arcing  is  in  connection 
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with  the  gas  recirculation  in  the  laser,  resulting  in  the  buildup  of  the  attaching  oxygen  compo¬ 
nent  in  the  laser  gas  mixture.  These  observations  could  also  be  explained  within  the  two-step 
ionization  model.  Although  the  details  of  the  kinetics  in  a  high  pressure  molecular  discharge  are 
complex,  the  concept  of  two-step  ionization  as  the  dominant  ionization  source  in  the  discharge  is 
plausible.  The  basic  reason  stems  from  the  fact  that  the  electron  temperature  in  these  discharges 
is  relatively  low  (typically  1-2  electron  volts)  and  very  few  electrons  in  the  tail  of  the  distribu¬ 
tion  function  can  directly  ionize  the  ground  state,  while  the  ionization  of  a  variety  of  available 
excited  and  metastable  states  requires  much  lower  energies  and  is  accessible  by  the  bulk  of  the 
electron  energy  distribution.  It  is  precisely  this  reason  which  makes  the  further  development  of 
the  multi-step  ionization  discharge  theory  necessary. 

2.2  The  Steady  State  Solution 

Consider  a  system,  shown  schematically  in  Fig.  2.1  with  a  DC  voltage  source  V„.  The 
voltage  is  distributed  between  the  external  resistor  R  and  the  plasma,  so  that 


V0  =  RI  +  LE 


(2.1) 


where  I  is  the  discharge  current,  L  is  the  plasma  length  and  E  is  the  electric  field  in  the  plasma, 
which  is  assumed  to  be  uniform  at  all  times.  In  writing  Eq.(2.1)  it  is  assumed  that  the  cathode 
fall  voltage  is  negligible  compared  to  LE  (a  situation  typical  of  high  ■pressure  discharges). 
Three  regimes  of  operation  will  be  investigated: 

A.  Stiff  Voltage  Case  -  Here  R  ~  0  and  thus  the  field  in  the  discharge 

V0 

E  =  — -  =  Const  (2.2) 

L* 


is  not  affected  by  a  possible  time  variation  of  discharge  characteristics,  such  as  the 
electron  density,  temperature  etc. 

B.  Stiff  Current  Case  -  In  this  case  R  — ►  oo  and  V0  — *  oo,  so  that  the  current 


/  ~ 


Vo 

—  =  const 
R 


(2.3) 


is  time  independent  and  fully  described  by  the  external  circuit 
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C.  Finite  External  Impedance  Case  -  In  this  intermediate  case  (0  <  R  <  oo)  both  the 


current 


m  = 


Vo  -  LE(t) 
R 


(2.4) 


and  the  electric  field  E(t)  depend  on  the  details  of  the  time  evolution  of  the  discharge. 
At  this  point  it  is  necessary  to  discuss  the  kinetics  of  the  self-sustained  discharge.  The 
densities  of  various  species  in  the  discharge  will  be  assumed  to  be  spatially  uniform,  and  the 
discharge  is  described  by  the  following  system  of  rate  equations 


dne 

dt 

dnm 

dt 


=  vnmnt  -  anl  -  (3nt 

=  ‘ynerta - 

r 


(2.5) 

(2.6) 


Here  nt  is  the  electron  density,  governed  by  the  ionization  of  an  excited  or  metastable  state  having 
a  density  nm  and  by  recombination  and  attachment  losses  (v,  a  and  /?  are  the  corresponding 
rate  constants).  The  density  of  the  excited  or  metastable  state  is  defined  by  the  balance  between 
the  processes  of  the  excitation  by  electron  impact  of  atoms  from  the  ground  state  (of  density 
na)  and  the  collisional  quenching  by  the  neutral  gas  (with  characteristic  time  constant  r).  The 
time  scales  r,  (ane)  and  /3-1  are  assumed  to  be  much  shorter  than  the  acoustic  times  and  thus 
the  neutral  gas  temperature  and  density  are  constants. 

Equations  (2.5)  and  (2.6)  yield  the  steady  state  in  which  all  the  discharge  characteristics  are 
constant  and 


r?m  r?m0  —  T a  (2.7) 

=  nt0  =  - - — - -  (2.8) 

(i'oT »Tna  -  a0) 

The  question  of  whether  this  steady  state  can  be  reached  for  a  given  set  of  parameters  of  the 
external  circuit  (V0,  R  and  L)  will  be  considered  next.  The  answer  to  this  question  depends  on 
the  regime  of  operation: 

A.  Stiff  Voltage  Case  -  Here  E  =  E0  is  given,  which  defines  all  the  rates  in  Eq.  (A. 8)  so 
that  neo  is  known.  This,  in  turn,  defines  the  current  density  j  =  ent0Vi0  in  the  plasma 
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( Vdo  being  the  drift  velocity  of  the  electrons).  Thus  the  steady  state  exists,  provided  of 
course  (see  Eq.  (2.8))  the  quantity 

Q(E0)  =  VoioTria  -  Q0  (2.9) 


is  positive.  This  condition  restricts  the  accessibility  of  the  steady  state  and  needs  a 
further  discussion.  Consider  the  dependence  of  Q  on  E0.  Typically,  the  excitation  and 
ionization  rates  (y0,  v0)  are  increasing  functions  of  the  electric  field,  while  the  attachment 
rate  a0  normally  decreases  with  E0.  Furthermore,  J°ZoQ  —  ~a(E0~ o)  <  0  and  with 
the  increase  of  E0,  Q  — >  v0~i 0Tna  >  0-  Therefore  there  exists  a  single  value  of  the 
electric  field  E0  =  E *  at  which  Q(E*)  =  0.  The  steady  state  operation  of  the  discharge 
under  the  stiff  voltage  conditions  ( E0  is  given)  is  thus  only  possible  when 


Eq  ~ 


E * 


(2.10) 


B.  Stiff  Current  Case  -  Here  the  current  I0  is  defined  by  the  external  circuit.  On  the  other 
hand  from  Eq.  (2.8) 

I0  =  evdori'0S  =  =  E(E0)  (2. 1 1) 

Q(E0 ) 


where  5  is  the  cross  section  area  of  the  discharge.  Equation  (2.1 1)  now  defines  the  steady 
state  electric  field  value  E0.  A  graphical  illustration  of  the  solution  of  this  equation  is 
given  in  Fig.  2.2,  where  the  qualitative  dependence  of  function  F  on  E  is  shown  based 
on  our  earlier  discussion  of  function  Q(E).  A  possible  minimum  of  F  at  E  =  Em  is 
also  shown.  Such  a  minimum  may  exist  at  high  values  of  E0  at  which  the  ionization 
and  excitation  rates  grow  slower  than  the  drift  velocity  vd  as  the  electric  field  increases. 
It  can  be  seen  in  Fig.  2.2  that  the  steady  state  exists  for  any  current  I0,  provided  it  is 
larger  than  Im.  It  is  noteworthy  that  the  operation  near  E'0  requires  rather  high  currents. 

C.  Finite  Impedance  Case  *  In  this  intermediate  case  Eq.(2.11)  is  replaced  by  (see  Eq. 
(2.4)  for  the  current) 


-  LE0 


=  F(E0) 


(2.12) 


R 
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The  graphical  solution  of  this  equation  for  E0  is  shown  in  Fig.  2.3.  One  observes  that 
for  given  values  of  V0  there  exists  a  critical  impedance  Rcr  such  that  when  R  <  Rcr 
two  steady  state  solutions  {E0\  and  £o2)  exist  for  the  electric  field,  while  for  R  >  Rcr 
no  steady  state  exists  at  all.  For  R  =  Rcr  there  is  a  single  steady  state.  The  value  of 
Rct  can  be  found  from 

L  =  —RCTF'(Ecr)  (2.13) 


where  F'  =  dF/dE0  and  Ecr  is  given  implicitly  by 


Vo  _  F(Ecr) 
L  F'(Ecr) 


(2.14) 


This  completes  the  analysis  of  the  possibility  of  a  steady  state  solution  at  various  external 
circuit  conditions. 

2.3  Numerical  Example 

At  this  stage  the  stability  of  the  discharge  at  various  operating  conditions  will  be  analysed. 
The  analysis  will  be  started  by  presenting  the  dependence  of  function  F(E)/S  on  the  electric 
field  (see  Eq.  (2.11)).  This  dependence  is  shown  in  Fig.  2.4  for  the  case  of  interest.  One 
observes  that  E *  ~  8.05  kV/cm,  so  that  the  steady  state  discharge  operation  is  possible  for 
8.05  kV/cm  <  E  <  15  kV/cm  (see  the  discussion  in  Section  2.2  and  Figs.  2.2  and  2.3).  The 
graphical  solution  for  the  electric  field  E0  in  the  steady  state  discharge  is  also  shown  under  the 
following  conditions:  V0  =  9  kV,  RS  =  400  ftcm2,  L  =  1  cm,  r  =  1  /. /sec  and  (30  =  105  sec-1. 
The  steady  state  electric  field  is  8.45  kV/cm  which,  by  using  Eq.  (2.8),  and  rate  coefficients 
appropriate  for  the  1  atm,  3:2:1  (He:N2:CG2)  gas  mixture  (determined  previously  as  function  of 
E),  yields  the  steady  state  plasma  density  nco  =  1.07  •  1012cm-3.  Prior  to  studying  the  stability 
of  the  discharge  at  these  conditions,  the  stiff  voltage  situation  is  considered.  Previously  reported 
analysis  showed  that  the  discharge  is  unstable  in  the  stiff  voltage  case  for  jd0  ^  0.  This  unstable 
behavior  is  demonstrated  in  Figs.  2.5  and  2.6,  where  the  numerical  solutions  of  Eqs.  (2.5)  and 
(2.6)  are  shown  for  E0  =  8.45  kV/cm  =  Const(f)  (stiff  voltage  condition)  and  two  different  initial 
values  for  the  electron  density.  The  electric  field  E0  defines  various  rates  in  Eqs.  (2.5)  and 
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Figure  2.3:  Graphical  solution  for  the  electric  field  Eo  in  the  discharge  under  finite  impedance 
conditions.  The  straight  lines  describe  the  function  (Vo  -  LEo)/R  for  different 
values  of  the  external  impedance  R. 
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TIME  (^sec) 


Figure  2.5:  The  time  evolution  of  the  electron  and  metastable  densities  in  a  stiff  voltage 

discharge.  The  initial  electron  density  (n^O)  =  5xlOn  cm3)  is  lower  than  the  steady 
state  density  (neo  -  1.07X1012  cm3)  at  the  given  voltage  (Vo  =  8.45  kV). 
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TIME  (jisec) 


Figure  2.6:  The  time  evolution  of  the  electron  and  metastable  densities  in  a  stiff  voltage 

discharge  with  attachment.  Initial  electron  density  (2x  1012  cm  3)  is  higher  than  the 
steady  state  density  (no,  =  1.07X1012  cm  3)  at  the  given  voltage  (Vo  =  8.45  kV). 
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(2.6)  which  remain  constants  in  the  calculations,  consistent  with  the  stiff  voltage  conditions. 
Figure  2.5  shows  the  case  when  ne(0)  =  nm( 0)  -  5  ■  10ncm-3.  This  initial  density  is  lower 
than  the  above  mentioned  steady  state  density  (ne0  =  107  •  1012cm-3).  Figure  2.5  shows  that 
the  electron  density  decays  in  time,  in  this  case,  and  the  discharge  is  quenched.  Figure  2.6,  in 
contrast,  shows  the  case  of  ne(o)  =  nm(o)  =  2  •  1012  cm-3  which  is  larger  than  the  steady  state 
election  density  r?.,..  nc  then  grows  indefinitely  in  time,  verifying  the  instability.  One  finds  that 
Figs.  2.5  and  2.6  represent  the  general  time  evolution  of  the  discharge  in  the  stiff  voltage  case, 
depending  on  whether  ne(o)  is  larger  or  smaller  than  neo. 

At  this  point  it  is  appropriate  to  proceed  to  the  numerical  simulation  of  the  finite  external 
impedance  discharge  case,  again  setting  V0  =  9  kV/cm,  L  =  1  cm,  r  =  1  //sec  and  /?  =  105  sec-1 
=  Const(t).  Nevertheless,  in  contrast  to  the  stiff  voltage  case,  the  rate  Eqs.  (2.5)  and  (2.6)  have 
now  been  solved  simultaneously  with  the  circuit  Eq.  (2.14).  The  code  also  makes  use  of  the 
dependencies  of  the  C02  gas  mixture  rate  constants  determined  previously.  The  predictions  of 
the  code  are  plotted  in  Figs.  2.7,  2.8  and  2.9.  Figure  2.7  shows  the  temporal  variation  of  ne  and 
E  for  the  circuit  impedance  of  RS  =  400  ficm2.  The  initial  electron  density  was  assumed  to  be 
1012cm-3.  One  observes  that  the  electron  density  increases  by  about  10%  to  1.07  •  1012cm-3 
and  the  electric  field  decreases  from  its  initial  value  of  9  kV/cm  to  8.45  kV/cm.  These  are  the 
steady  state  electron  density  and  the  electric  field  under  the  above  conditions.  Figure  2.8  shows 
the  effect  of  lowering  the  external  impedance  to  RS  =  1  cm2.  Under  these  conditions,  as  can 
be  seen  in  the  figure,  the  discharge  has  stiff  voltage  characteristics  during  the  first  ~  10  //sec 
of  the  discharge.  It  is  unstable  as  predicted  by  the  analysis  and  the  electron  density  suddenly 
increases  by  almost  three  orders  of  magnitude  after  10  //sec.  The  discharge  becomes  stable  at 
this  large  value  of  ne  ~  8  •  10Hcm-3.  However,  in  practice,  any  initial  spatial  nonuniformities 
in  ne  amplified  in  time  by  the  instability  could  result  in  streamer  formation  during  the  fast 
transition  period.  We  shall  consider  the  details  of  such  spatial  instabilities  in  the  next  chapter  of 
this  report  Finally,  Fig.  2.9  illustrates  a  typically  situation  in  a  discharge  with  an  above-critical 
external  impedance  ( R  >  RCT,  see  the  discussion  in  the  previous  section).  In  the  case  considered 
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TIME  (^sec) 

Figure  2.7:  Curves  of  ne  and  E  for  a  volumetrically  stable  discharge.  The  parameters 

are:  ne(0)  =  1012  cm  3,  Vo  =  9  kV/cm,  RS  =  400  Q  cm2, 1  =  1  cm,  t  =  1  nsec, 
P  =  105  sec1. 
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Figure  2.8:  Curves  of  ne,  and  E  for  a  low  external  impedance  discharge  (RS  =  1  W  cm2). 

The  parameters  are  ne(0)  =  10^  cm-3,  Vq  =  9  kV/cm,  L  =  1cm,  t  =  1msec  and  b  =  10^ 
sec-1.  The  discharge  is  unstable  during  the  initial  evolution  phase  (t  <  15  msec). 
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Figure  2.9:  The  curves  of  ne  and  E  for  a  discharge  with  an  above-  critical  impedance 

(RS  =  800  ft  cm2).  All  other  parameters  are  the  same  as  in  Figs.  3.10  and  3.1 1.  No  steady 
state  exists  for  RS  >  R^S  =  520  ft  cm2  and  the  discharge  is  quenched. 
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(see  Fig.  2.4)  Rcr  ~  520  ficm2,  while  the  results  presented  in  Fig.  2.9  correspond  to  the  case 

R  =  800  Qcm2.  No  steady  state  exists  in  this  case,  as  predicted  by  the  theory,  and  one  observes 

in  Fig.  2.9  that  the  discharge  is  quenched. 

2.4  Stability  Conclusions 

(i)  We  have  considered  the  volumetric  stability  dynamics  of  discharges  dominated  by  two  step 
ionization.  In  the  presence  of  an  attacher  these  discharges  are  shown  to  be  potentially  unsta¬ 
ble  on  much  shorter  time  scales  than  those  characteristic  of  conventional  thermal  instabilities. 

(ii)  The  theoretical  modeling  of  the  discharge  included  the  discharge  kinetics  coupled  to  the 
external  electric  circuit  and  its  response  to  varying  current  conditions.  A  number  of  dis¬ 
charge  configurations,  such  as  stiff  voltage,  stiff  current  and  finite  external  impedance  were 
considered  in  searching  for  a  steady  state  stable  regime  of  operation. 

(iii)  In  the  stiff  voltage  discharge  [the  electric  field  in  the  discharge  E  =  const  (t)]  the  steady 
state  exists  for  E  >  E\  where  E *  is  the  value  at  which  the  characteristic  function  Q(E) 
[see  Eq.  (2.9)]  vanishes.  Nevertheless,  even  if  the  steady  state  is  reached,  the  discharge 
is  volumetrically  unstable  under  the  stiff  voltage  condition.  Practically,  if  at  t  =  0  the 
initial  electron  density  n<( 0)  >  neo  ( nc0  being  the  steady  state  density)  the  current  density  j 
increases  indefinitely  on  a  finite  time  scale,  while  for  ne( 0)  <  nto  the  discharge  is  quenched 
(j(+)  — >  0).  The  linear  instability  temporal  rate  is  typically  proportional  to  the  attachment 
rate  0  (if  0  is  sufficiently  small). 

(iv)  The  stiff  current  discharge  [j  =  const(t)]  has  a  steady  state  at  all  currents  and  the  discharge  is 
volumetrically  stable,  provided  the  characteristic  function  F(E)  [see  Eq.  (2.1 1)]  is  decreasing 
function  of  the  electric  field. 

(v)  In  the  finite  external  impedance  discharge  case  we  found  the  existence  of  a  critical  impedance 
value  Rct  such  that  there  exist  two  steady  state  when  0  <  R  <  Rcr  and  no  steady  state  is 
possible  when  R  >  Rcr.  Furthermore,  only  if  Rcr  is  large  enough,  the  above  mentioned 
steady  state  solutions  are  stable  in  a  part  of  the  interval  [0,  Rcr ]. 

(vi)  The  theoretical  predictions  were  illustrated  in  the  example  of  an  atmospheric  C02  laser 
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discharge  model.  The  numerical  solutions  of  the  kinetic  and  external  circuit  equations, 
underlying  the  problem  in  this  case  were  found  in  an  agreement  with  the  theory. 

2.5  Numerical  Simulations  of  The  Two  Dimensional  Code 

SRL  has  developed  a  code  based  on  the  two  dimensional  RF  discharge  theory  presented  in 
a  previous  report.  The  basic  equation  in  rectilinear  coordinates  is 


d2E  i4ne2neo  LO 

dx 2  me2  vcou 


(2.15) 


subject  to  the  boundary  conditions 


dE 

dx  la 


=  0  and 


dE 

dx 


1 1—  w 


2 ujL 


where  I0  is  the  total  current  per  unit  length.  Equation  (2.15)  is  valid  for  a  slab  geometry  i.e. 
where  the  2  dimension  is  much  longer  than  x-  dimension.  Such  a  geometry  conforms  more  to  the 
C02  lasers  than  the  radial  geometry.  It  is  assumed  that  the  frequency  of  oscillation  /  >>  cme, 
hence  the  electron  density  will  remain  essentially  constant  during  a  period  of  oscillation.  This 
allows  one  to  use  time  averaged  rates  for  a,  v  and  (ov).  These  time  averaged  rates  are  defined 
by  the  following  integral 

*7^  r 

l  =  ^  J  £(|  E  |  sinwOdf  =  ^  J  £(E0sins)</s  =  |(|  E  |)  (2.16) 

o  0 

where  £  =  v,  a  and  (av). 

In  Eq.  (2.16)  |  £  |  is  just  the  amplitude  of  the  electric  field 


|  E  |=  {(ReE)2  +  (ImE)2}i/2 


(2.17) 


For  efficient  simulation  one  has  to  obtain  analytic  expressions  for  v,  a  and  (a v)  as  functions 
of  |  E  |.  It  is  also  convenient  to  introduce  a  dimensionless  coordinate  p  =  x/w  in  which  case 
Eq.  (2.15)  can  be  rewritten  as 

d2  (  w \ 2 

wE+2i(j)  E=°  <2I8> 
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where  the  skin  depth  6  is  given  by 


S  -  (S—\ 1/2  =  (  cmuc°“  ^ 
\2naujJ  \2n  e2ntu  J 


1/2 


It  is  assumed  that  the  discharge  width  w  «  6  and  so  one  can  solve  Eq.  (2.15)  by  a 
perturbation  expansion  by  introducing  the  small  parameter  A  =  ( w/6 )2.  The  electric  field  can 
be  expanded  as 

E  =  E0  +  E\  +  Ei . 


where  the  terms  are  ordered  in  powers  of  A  then  Eq.  (2.15)  gives  in  various  orders 

d2Ea 


Op 2 


=  0 


d2E1 
dp 2 


+  2iAE0  =  0 


. etc. 


(2.19) 

(2.20) 
(2.21) 


Equation  (2.19)  predicts  a  constant  value  for  E0  which  after  substitution  into  (2.20)  results  in 


p  p 

Ei(p,t)  =  —2iEo  J  dp'  J  A (p")dp" 


(2.24) 


Similarly 


E2(p,t)  =  -2 i  j  dp'  J  dpE\Adp” 

O 

etc.  The  SRL  code  solves  the  first  order  equation  given  by 

E(p,  t)  =  E0  ^1  -  2 «■  jT  dp'  J  A(p",  t)dp" 

where  A  is  found  by  solving  the  following  rate  equations: 

dnt 


(2.25) 


(2.26) 


dt 


=  vnmnt  —  an*  —  0ne  +  DaV*ne 

dnm  nm  2 

— -  =7 nt - +  Dmv  nm 

dt  t 
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Da  and  Dm  are  the  ambipolar  and  metastable  atomic  diffusion  coefficients  respectively.  Dm  is 
of  the  order  l/Zv^X  ~  102  cm2/sec/Torr  while  Da  may  probably  be  100  times  larger  ( Da  ~  104 
cm2/sec  Torr).  The  constant  E0  can  be  found  from  the  boundary  condition  on  the  total  current 


ae 

dp 


\p-\ 


2iu>In 


=  -2  iE, 


J  A  {p\t)p'dp' 


or 

E0  =  ^/  f\(p\t)p'dp'  (2.29) 

The  numerical  code  results  will  be  presented  next.  The  particular  slab  geometry  analysed  is 
shown  in  Fig.  2.10.  The  discharge  is  assumed  to  be  2.4  cm  wide  in  the  x-  direction  and  long 
in  the  z  direction.  The  electric  field  is  applied  in  the  y  direction.  The  origin  of  the  abscissa  is 
located  at  the  center  of  the  discharge  i.e.  the  discharge  extends  between  -1.2  cm  to  1.2  cm  in 
the  x-  direction.  The  stiff  current  drive  will  be  discussed  first  since  it  was  the  most  stable. 

2.5.1  A  Stiff  Current  Case 

Figure  2.1 1  shows  the  temporal  and  spatial  variation  of  ne  for  the  discharge  between  0  <  x  < 
1.2  cm.  The  discharge  is  assumed  to  be  symmetric  hence  -  1.2  cm  <  x  <  0  will  just  be  the 
mirror  image  about  the  ordinate  of  the  plots  shown  in  Fig.  2.11.  Also  it  was  assumed  that 
the  total  current  was  2.4  A/cm  and  the  applied  electric  field  was  400  V/cm  and  the  C02  laser 
mixture  was  3/2/1  at  a  total  pressure  of  30  Torr.  The  initial  spatial  electron  density  was  assumed 
to  be  slightly  greater  than  10n  cm-3  with  a  non  uniformity  of  5%  between  the  centre  of  the 
discharge  and  the  edge.  The  attachment  rate  for  electrons  was  assumed  to  be  50  // s.  The  spatial 
variation  of  the  electrons  is  shown  every  5  ps.  The  initial  spatial  profile  is  labeled  1.  The  curve 
labeled  2  is  spatial  profile  5  //s  later  etc.  The  electron  density  increased  by  20%  in  the  first  10 
//s  and  then  decrease  by  10%  during  the  next  50  ps.  Also  the  non  uniformity  at  the  end  of  45 
ps  increased  to  about  10%.  Figure  2.12  shows  the  results  of  the  next  45  ps.  The  curves  in 
Figure  2.12  are  spaced  every  10  //sec.  The  curve  labeled  1  in  Fig  2.12  is  the  same  as  the  curve 
labeled  10  in  Fig.  2.11.  The  curve  labeled  2  shows  the  spatial  variation  of  ne  10//S  later,  etc. 
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Figure  2  10:  Schematic  of  discharge  geometry 
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ELECTRON  DENSITY  AT  DIFFERENT  TIMES 
(At  =  5  }J.S) 


Figure  2.11:  Plot  showing  the  spatial  and  temporal  variation 
of  the  electron  density  ne  .  Curve  1  is  the  initial 
spatial  variation  of  ne  at  t  =  0.  Curve  2  is  the 
spatial  variation  of  ne  at  t  =  5  (is,  etc.  The  power 
supply  is  assumed  to  be  a  stiff  current  source. 
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ELECTRON  DENSITY  AT  DIFFERENT  TIMES 


Figure  2.12:  The  curve  labled  1  is  the  same  as  the  curve  labeled  10 
in  Figure  2.1 1.  This  figure  shows  the  spatial  variation 
of  ne  every  10  ps. 
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The  dominant  effect  is  the  continual  increase  in  the  non  uniformity.  The  uniformity  at  the  end 
of  85  ^s  is  25%. 

The  result  shown  in  Fig.  2.11  and  2.12  show  for  the  conditions  assumed  the  discharge  was 
stable  for  about  100  us.  However  the  spatial  non- uniformity  increases  temporally  indicating  the 
discharge  will  eventually  be  unstable  and  arcing  will  occur.  From  the  many  code  runs  analysed  it 
is  clear  that  a  spatial  non-uniformity  and  a  finite  attachment  rate  results  in  an  unstable  discharge 
even  for  a  stiff  current  source  and  the  key  for  stable  operation  is  to  obtain  as  uniform  an  initial 
electron  density  as  possible  and  to  keep  the  electron  attachment  rate  as  small  as  possible. 

2.5.2  Stiff  Voltage  Source 

Figures  2.13  and  2.14  show  the  results  of  a  stiff  voltage  source  for  the  same  conditions 
discussed  previously.  In  Fig.  2.13  the  initial  electron  density  was  again  chosen  to  be  about 
1.1  xlO11  cm'3.  The  electron  density  increases  by  over  an  order  of  magnitude  in  30  ^/s.  In  this 
figure  the  spatial  variation  in  the  electron  density  is  shown  every  3  /rs  as  opposed  to  every  5 
US  in  Fig.  2.11  and  every  10  fi s  in  Fig.  2.12.  This  discharge  is  clearly  unstable  and  it  will 
result  in  an  arc  well  before  the  100  /ts  pulse  length.  Figure  2.14  shows  the  results  of  initializing 
the  electron  density  to  1010  cm'3.  In  this  case  the  discharge  is  clearly  quenched.  These  results 
are  reminiscent  of  the  one  dimensional  stiff  voltage  results  shown  in  Figs.  2.8  and  2.9  where 
discharge  is  unstable  and  it  proceeds  to  an  arc  or  its  is  quenched  depending  on  the  initial  electron 
density.  So  once  again  it  is  clear  that  the  stiff  current  source  is  more  stable  than  a  stiff  voltage 
source. 
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TIME 

TEMPORAL  SEPARATION 
BETWEEN  CURVES 
=  3  usees 


Figure  2.13:  Same  as  Figures  2.11  and  2.12  except  that  power  supply 
is  assumed  to  be  a  stiff  voltage  source  and  the  spatial 
variation  of  ne  is  shown  every  3  its. 


2-24 


SCIENCE  RESEARCH  LABORATORY 


1.1e+10 


I.O0+IO 


9.06+9 


8.0e+9 


7.00+9 


6.O0+9 


5.00+9 


4.0e+9 


3.00+9 


|  TIME 

TEMPORAL  SEPARATION 
BETWEEN  CURVES 
=  3  psecs 


0.0  0.2  0.4  0.6  0.8  1.0  1.2 

DISTANCE  FROM  CENTER  (cm) 


Figure  2.14:  Same  as  Figtire  2.13  except  that  ne 
was  initialized  at  10 10  cm'3 
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SECTION  3 


3.0  STEADY-STATE  CONDUCTION  COOLING 


In  a  repetitive  pulsed  C02  laser,  the  lasing  gas  must  be  cooled  to  prevent  thermal  population 
of  the  lower  laser  level  and  a  decrease  in  laser  efficiency.  At  a  repetition  rate  of  v  pps,  an 
extraction  of  W  (J/liter),  a  volume  of  V  (liters)  and  an  electrical  efficiency  of  r?,  the  amount  of 
heat  Q  which  must  be  carried  away  from  the  laser  to  maintain  it  at  constant  temperature  is 


Q  = 


uWV 


(3.1) 


For  example,  for  v  =  50  pps,  W  =  1.67  J/liter,  V  =  6  liters  and  r?  =  0. 167,  values  typical  of  a  C02 
laser  operating  at  a  pressure  of  30  Torr,  one  obtains  Q  «  3.0  kW.  Convecting  away  this  much 
heat  typically  requires  a  flow  loop,  which  increases  system  weight  substantially.  In  this  section 
the  alternative  approach  of  conduction  cooling  will  be  addressed  instead  as  an  advanced  concept 
to  minimize  the  total  weight  of  the  laser  head.  As  shown  in  Fig.  3.1,  the  heat  may  be  conducted 
away  by  the  sidewalls  along  a  direction  perpendicular  to  the  discharge  direction.  Conduction 
cooling  is  a  viable  option  if  the  laser  is  operated  at  a  low  pressure  and  a  correspondingly  low 
extraction  energy  density.  Lower  extraction  energy  density  and  increased  laser-head  volume 
result  in  a  lower  overall  laser  head  weight  when  they  are  linked  to  the  significant  weight  and 
volume  reductions  made  possible  by  the  elimination  of  the  flow  loop  and  associated  blowers, 
acoustic  damping  system  and  heat  exchangers. 

3.1  Conduction  Equations 

For  a  repetitively  pulsed  C02  laser  device  the  temperature,  density  and  pressure  of  the 
lasing  gas  all  vary  during  each  pulse.  A  repetitively  pulsed  C02  laser  applicable  to  laser  radar 
is  expected  to  have  pulses  of  width  ~  30  //sec.  These  pulses  are  so  short  that  conduction  within 
the  laser  gas  may  be  ignored  during  each  pulse.  The  equations  coupling  the  temperature,  density 
and  pressure  may  then  be  reduced  to  a  third-order  equation  for  the  density  disturbance  driven  by 
the  pulsed  heating.  The  finite  acoustic  speed  limits  the  response  of  density  disturbances  of  scale 
sizes  greater  than  ~  1  cm  during  each  pulse.  Between  pulses  the  situation  is  quite  different, 
however.  The  pulses  are  far  enough  (~  20  msec)  apart  for  gaseous  conduction  of  the  heat  to 
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Figure  3.1:  Schematic  of  conduction-cooled  CO  laser  head 
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be  significant.  Since  the  rate  of  heat  conduction  is  proportional  to  the  temperature  gradient, 
small-scale  temperature  gradients  will  tend  to  conduct  away,  leaving  behind  predominantly  the 
largest-scale  gradients.  After  a  large  number  of  such  pulses,  the  laser  will  reach  a  quasi-  steady 
state,  in  which  the  temperature  and  density  are  periodic  in  time  with  the  period  of  the  electrical 
pumping.  In  this  limit,  the  effect  of  each  individual  electrical  pulse  on  the  beam  quality  is  not 
expected  to  be  severe.  However,  the  cumulative  effect  of  a  large  number  of  pulses,  balanced 
against  the  conduction  loss,  is  expected  to  lead  to  steady-state  density  and  temperature  profiles 
which  significantly  influence  the  wavefront  of  the  laser  beam.  This  steady-state  effect  may  be 
calculated  approximately  by  ignoring  the  pulsed  nature  of  the  excitation  and  treating  it  as  if  the 
heating  is  uniform  in  time,  with  the  same  average  rate  as  that  of  the  actual  repetitively  pulsed 
heating. 

For  heating  which  is  uniform  in  time,  the  equation  for  the  temperature  T,  including  the 
effect  of  thermal  conduction  of  the  gas,  is 


-“V’T=S 


(3.2) 


where  cp  is  the  specific  heat  (at  constant  pressure)  per  unit  volume,  k  is  the  thermal  conductivity 
of  the  gas,  and  S  is  the  heating  rate  per  unit  volume.  Eq.  (3.2)  may  be  written  in  the  alternative 
form  1 

^-£>V2T  =  -  (3.3) 

at  cp 

where  D  =  k/cp  is  the  heat  diffusivity,  which  scales  as  the  product  of  the  acoustic  speed  and  the 
mean-free  path  between  collisions.  It  is  thus  inversely  proportional  to  the  gas  density.  Therefore, 
if  the  heating  source  term  S/cp  is  kept  constant  as  the  density  is  decreased  from  1  atm,  the  effect 
of  heat  conduction  increases  significantly.  This  may  also  be  seen  directly  from  Eq.  (3.2):  if  one 
decreases  the  density  and  simultaneously  decreases  S  proportionally,  then  since  k  is  independent 
of  density  the  conduction  term  takes  on  increasing  significance  as  the  density  is  lowered. 

An  estimate  of  the  efficacy  of  conduction  cooling  for  relevant  laser  conditions  may  be  made 
directly  from  Eq.  (3.2).  Consider  a  laser  channel  of  dimensions  2.4  x  20  cm2,  with  the  surfaces 
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separated  by  d  =  2.4  cm  maintained  at  constant  temperature  T0.  For  a  3:2:1  He:N2:C02  mix 
at  400°K,  we  have  k  =  1  x  10-3W/cm°K.  In  this  geometry  the  cooling  is  dominated  by  the 
effect  of  the  walls  spaced  2.4  cm  apart.  Neglecting  the  other  walls,  assuming  a  spatially  uniform 
heat  deposition,  and  taking  the  steady-state  limit  (T  independent  of  time),  one  may  integrate  Eq. 
(3.2)  analytically  to  obtain 

Wv 

A T(x)  =  —x(d  -  x)  (3.4) 

2  KTJ 

where  x  is  the  distance  from  one  of  the  walls.  At  v  =  50  pps,  W  =  1.67  x  10“3  J/cm3,  r;  =  0.167, 
Eq.  (3.4)  predicts  a  steady-state  temperature  rise  of  ~  360  K  halfway  between  the  walls. 

The  1.67  J/liter  extraction,  16.7%  efficiency  and  50  pps  repetition  rate  will  be  taken  to  be  the 
baseline  parameters  of  a  C02  laser  amplifier  to  be  used  for  laser  radar.  A  preliminary  analysis 
indicates  that  the  d  =  2.4  cm  spacing  gives  adequate  conduction  cooling  at  v  =  50  pps;  therefore, 
unless  otherwise  noted,  this  spacing  will  be  tentatively  considered  a  baseline  parameter,  as  well. 
At  lower  repetition  rates,  a  larger  wall  spacing  can  be  tolerated. 

3.2  Thermal  Effect  of  Pulsed  Nature  of  Heating 

One  can  also  solve  Eq.  (3.2)  for  the  temperature  as  a  function  of  time,  assuming  that  the 
energy  deposition  in  the  medium  is  uniform,  and  that  the  changes  take  place  at  constant  pressure. 
One  obtains  a  maximum  temperature  rise  A Tm  given  by2 


A  Tm  =  A  Tp 


(e®  —  1  3(e9,>  —  1)  ) 


(3.5) 


In  Eq.  (3.5)  A Tp  =  W/cptj  is  the  temperature  rise  due  to  energy  deposition  by  a  single  pulse 
and  the  dimensionless  constant  v  is  given  by 


7T\  2 


V  = 


VC. 


G) 


(3.6) 


where  cp  is  the  specific  heat  at  constant  pressure.  For  d  =  2.4  cm,  one  obtains  v  «  0.83, 
A Tp  «  220  K,  and  ATm  «  420  K. 

Finally,  it  is  useful  to  determine  the  rate  at  which  temperature  nonuniformities  of  differ¬ 


ent  scale  sizes  are  dissipated  by  thermal  conduction  between  pulses.  Let  ^2o(r)  represent  the 
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4 


nonuniform  temperature  rise  at  the  end  of  a  pulse.  Expanding  the  temperature  in  Eq.  (3.2)  in 
spatial  Fourier  components  yields  during  the  time  between  pulses 


dmk)+Kk26^  0 

at 

with  the  solution 

6T(k)  =  <5t0a-)e~Kp</Cp 


At  a  repetition  rate  of  v  Hz,  the  pulses  are  separated  by  a  time  t  =  \/u.  Therefore  the  residual 
temperature  rise  of  spatial  frequency  kn  =  nn/d  due  to  a  single  previous  pulse  is  given  by 


6Tn 

6Ton 


-  exp( 


— /C7T2n2 
cpd2  v 


) 


where  v  is  defined  in  Eq.  (3.6),  and  n  =  1,2...  is  an  integer  which  characterizes  the  spatial 
frequency  as  a  multiple  of  ix/d.  As  discussed  above,  v  «  0.83  when  d  =  2.4  cm.  Therefore,  the 
temperature  component  with  the  lowest  spatial  frequency  (n  =  1)  supported  by  the  channel  decays 
by  a  factor  of  ~  exp(0.83)  =  2.3  between  pulses.  The  component  with  the  next  lowest  spatial 
frequency  (n  =  2)  decays  by  a  factor  of  ~  exp(6)  =  28  between  pulses.  Thus  the  small-scale-size 
nonuniformities  quickly  dissipate  away  between  pulses,  and  to  an  excellent  approximation  only 
the  largest-scale  temperature  nonuniformity  needs  to  be  taken  into  account. 

3.3  Scaling  of  Conduction-Cooled  Laser 


Eq.  (3.4)  may  be  used  to  scale  the  laser  size  as  a  function  of  electrical  power  loading. 
The  steady-state  temperature  Tc  halfway  between  the  walls  may  be  expressed  in  terms  of  the 
electrical  power  loading  P'  =  Wvjrj  as 


Thus  the  temperature  rise  AT  is  proportional  to  the  power  loading  and  the  square  of  the  distance 
between  the  walls.  This  relationship  is  depicted  graphically  in  Fig.  3.2,  which  shows  the 
temperature  rise  as  a  function  of  the  power  loading  for  various  values  of  channel  width  (wall 
separation). 
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The  scaling  with  pulse  repetition  frequency  (PRF)  v  may  be  obtained  by  rewriting  Eq.  (3.7) 


in  the  form 


2  8k  Tc  —  To 


dl  = 


(3.8) 


QV/ri)  v 

In  Fig.  3.3  the  relationship  between  channel  width  and  PRF  is  plotted  for  various  values  of 
temperature  rise,  assuming  a  single-pulse  electrical  energy  loading  W/t]  of  10  J/lit. 

Table  3.1  lists  steady-state  temperatures  and  transient  temperature  rises  for  a  laser  head  with 
a  2.4-cm  channel  spacing  operated  at  various  pulse  repetition  rates  (PRF’s).  At  the  highest 
PRF  shown  (50  pps)  it  is  necessary  to  maintain  the  channel  walls  at  a  temperature  of  200°K 
(as  indicated  in  Table  3.1)  in  order  to  prevent  the  gas  from  reaching  temperatures  at  which  the 
lasing  efficiency  would  be  degraded.  The  200°K  temperature  is  still  well  above  the  boiling  point 
(~  160-170  K)  of  C02  over  the  estimated  pressure  range  (30-80  torr)  which  will  apply  during 
operation  of  the  laser. 


Table  I:  Temperature  Variation  for  2.4-cm-Channel  Laser  Head 


Steady  State  Temperature  (K) 

Transient  Temperature  Rise  (K) 

Repetition 
Rate  (Hz) 

At  Electrode 

At  Channel 
Center 

Peak 

Beginning  of 

Pulse 

10 

295 

370 

210 

0 

20 

295 

440 

250 

30 

50 

200 

560 

420 

220 

3.4  Impact  of  Conduction  Cooling  on  Beam  Quality  and  Spectral  Width 
in  the  Steady  State 

In  Section  3.1  it  was  shown  that  for  a  uniform  heat  deposition,  conduction  cooling  leads  to 
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PULSE  REPETITION  FREQUENCY  v  (pps) 


Figure  3.3:  Maximum  channel  width  versus  PRF  for 
various  steady-state  temperature  rises 


3-8 

SCIENCE  RESEARCH  LABORATORY 


a  steady  state  w<th  a  quadratic  temperature  variation  across  the  channel  [see  Eq.  (3.4)].  The 
quadratic  variation  of  the  temperature  will,  via  the  equation  of  state,  lead  to  a  corresponding 


variation  in  density  p: 


P  -  A  T0  + 


Wux(d  —  x y 


where  the  constant  A  may  be  determined  from  the  constraint  that  the  mean  value  of  p,  averaged 
over  x,  must  be  independent  of  the  heating.  The  main  aberration  coming  from  this  density 
variation  will  be  a  negative  cylindrical  focus  along  x,  i.e.,  perpendicular  to  the  channel  walls. 
The  associated  radius  of  curvature  R  may  be  calculated  from 


1  A  , 

—  =  — V2<3 
R  4tt  v 


(3.10) 


where  A  is  the  wavelength,  <f>  is  the  phase  aberration,  and  V2<£  is  evaluated  at  x  =  d/2,  halfway 


between  the  channel  walls.  The  relationship  between  phase  and  density  is 

,  2? rLl  /3g 

<t>  =  — — 

A  pA 


(3.11) 


where  jdg  is  the  Gladstone-Dale  constant  (the  refractive  index  minus  unity  at  standard  temperature 
and  pressure  [STP])  of  the  laser  gas,  pA  is  the  gas  density  at  STP,  and  L]  is  the  propagation 


length.  Combining  Eqs.  (3.9)-(3.11)  yields 

_L  ~  -h±R  p0  p 

R  ~  3  d?Pg pA  l+n 

where  p0  is  the  initial  density  of  the  gas,  and  p  is  given  by 

=  Tc  -  To  _  Wvd 2 
P  To  8  k  i]  Tq 


(3.12) 


(3.13) 


Here  T0  is  the  wall  temperature  and  Tc  is  the  temperature  at  the  center.  The  radius  of  curvature 
for  a  single  pass  through  the  amplifier  (Lx  =  125  cm)  when  po/pA  =  0.04  and  p  =  1  is 

R  «  120  m 

which  is  reasonably  large.  For  4  passes  through  the  amplifier,  with  no  correction  between  passes, 
the  radius  of  curvature  would  still  be  ~  30  m,  corresponding  to  ~  0.03  diopters.  The  cylindrical 
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focus  error  can  easily  be  compensated  statically  by  giving  the  laser  beam  incident  on  the  laser 
amplifier  a  cylindrical  focus  along  x.  If  the  amplifier  utilizes  an  expanding-wave  geometry  then 
the  optics  producing  the  expansion  (perpendicular  to  x)  will  be  in  place,  and  all  that  will  be 
needed  is  a  small  modification  of  these  optics  to  produce  the  relatively  small  focus  along  x.  Thus, 
the  cylindrical  focus,  which  is  the  dominant  aberration  due  to  the  steady-state  thermal  loading, 
can  easily  be  compensated  and  should  produce  no  appreciable  beam-quality  degradation. 

Higher-order  aberrations  due  to  the  thermal  loading  may  also  be  present  if  the  electrical 
heating  is  not  uniform  along  x.  These  may  be  estimated  by  taking  the  spatial  Fourier  transform 
of  Eq.  (3.2)  in  the  steady  state,  yielding 


T(k)  = 


S(k ) 
Kk 2 


(3.14) 


The  amplitude  of  the  component  of  the  temperature  of  spatial  frequency  k  is  thus  shown  to  be 
proportional  to  the  heating  nonuniformity  of  the  same  spatial  frequency  and  inversely  proportional 
to  the  square  of  the  spatial  frequency.  Therefore,  as  discussed  above,  small-scale  nonuniformities 
of  the  steady-state  heating  do  not  produce  strong  temperature  and  density  nonuniformities.  The 
reason  is  that  for  a  given  temperature  difference,  heat  is  conducted  more  rapidly  the  smaller  the 
scale  size  of  the  temperature  difference.  In  the  steady  state  the  dominant  density  disturbance  will 
always  be  the  low-order  cylindrical  focus,  which,  as  discussed  above,  is  easily  compensated.  If  it 
turns  out  to  be  necessary,  it  is  also  possible  to  compensate  some  of  the  higher-order  aberrations  by 
utilizing  a  static  adaptive  optic  which  may  be  adjusted  in  place  to  minimize  the  laser  divergence. 

It  should  be  noted  that  the  steady-state  temperature  and  density  profile  distributions  have  no 
effect  on  the  spectral  purity  of  the  laser.  Only  the  pulsation  of  the  heating,  which  causes  time- 
dependent  changes  in  the  phase,  can  broaden  the  spectrum  by  producing  chirp  and  frequency 
noise. 

A  preliminary  analysis,  sketched  above,  thus  indicates  that,  as  far  as  the  steady  state  is 
concerned,  conduction  cooling  to  channel  walls  is  a  viable  technique  for  carrying  away  the 
heat  generated  in  a  pulsed  C02  laser  amplifier  operated  at  a  low  pressure  (0.04  atm).  The 
heat  generated  by  each  electrical  pulse  is  sufficient  to  raise  the  temperature  of  the  laser  gas 
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significantly,  on  the  order  of  ~  300  °C.  However,  most  of  this  heat  conducts  away  to  the 
walls  between  pulses.  After  several  laser  pulses,  a  quasi-steady-state  temperature  and  density 
profile  appears.  It  is  comparable  with  the  profile  that  would  obtain  if  the  pulses  were  replaced 
by  continuous  heating  at  the  same  rate.  The  profile  is  dominated  by  a  quadratic  cylindrical 
foe  us  which,  for  four  passes  through  the  amplifier,  is  expected  to  be  ~  0.03  diopters.  This 
negative  focus  is  transverse  to  the  beam-expansion  direction,  and  is  easily  removed  with  a  slight 
modification  of  the  focussing  optics.  It  is  noteworthy  that  the  negative  focus  does  not  appear 
until  the  laser  amplifier  is  operated  for  several  pulses  so  that  it  reaches  its  steady-state  thermal 
balance.  Higher-order  aberrations  due  to  the  steady-state  heating  are  predicted  to  be  small, 
since  the  temperature  gradients  which  drive  them  are  strongly  dissipated  by  heat  conduction. 
However,  it  is  possible  to  remove  these  aberrations,  as  well,  by  utilizing  a  relatively  simple 
static  adaptive-optics  system.  The  steady-state  heating  should  have  no  effect  on  the  spectral 
purity  of  the  amplified  laser  beam. 
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SECTION  4 


4.0  ELECTRICALLY  GENERATED  DENSITY  DISTURBANCES 

The  electrical  energy  used  to  pump  the  laser  medium  will  unavoidably  have  some  degree  of 
spatial  nonuniformity.  This  will  result  in  the  generation  of  acoustical  and  entropy  disturbances, 
which  will  in  turn  produce  refractive-index  nonuniformities,  an  increased  spectral  width,  and 
beam-quality  degradation  in  the  laser.  Because  the  channel  walls  of  the  discharge-pumped  laser 
are  solid,  there  is  no  intrinsic  spatial  discontinuity  in  the  pump-energy  distribution,  such  as  that 
which  occurs  in  an  e-beam-pumped  laser.3  Still,  if  the  laser  medium  of  interest,  discharge-pumped 
by  ~  10  J/lit  in  a  ~  30-//sec  pulse,  were  used  as  an  oscillator,  it  would  be  extremely  difficult  to 
maintain  the  pumping  uniform  enough  to  make  the  oscillator  beam  quality  and  spectral  quality 
acceptable.  However,  as  will  be  shown,  the  beam-quality  degradation  and  frequency  noise 
are  not  so  severe  if  the  laser  medium  is  used  instead  as  an  amplifier.  It  is  expected  that  the 


dominant  source  of  the  discharge  nonuniformities  will  be  nonuniform  uv  preionization.  It  is 
thus  possible  to  make  the  laser  light  average  over  the  discharge  nonuniformity  by  ensuring  that 
the  uv-preionizing  structure  changes  along  the  direction  of  propagation  of  the  laser  light.  The 
distance  over  which  this  structure  must  change  significantly  in  order  to  ensure  good  beam  quality 
and  frequency  purity  can  be  determined  in  detail.  In  this  report,  a  preliminary  estimate  of  the 
requirements  for  good  beam  quality  will  be  made. 

4.1  Density  Response  to  Nonuniform  Pumping 

The  geometry  of  the  laser  is  shown  in  Fig.  4.1.  The  dominant  direction  of  discharge 
nonuniformity  is  expected  to  lie  along  the  direction  x  transverse  to  the  channel  surfaces.  Thus, 
the  contributions  of  nonuniform  pumping  to  density  variations  along  the  discharge  direction  y 


and  the  optic  axis  x  may  be  neglected,  and  to  a  good  approximation  the  perturbation  p  of  the 
density  satisfies  a  one-dimensional  wave  equation  of  the  form4 


d  (d2  2  d2  \ 

dt  l.3f2  “  3x2  J 


(4.1) 


where  Q(x,t)  is  the  heating  rate,  u  is  the  speed  of  sound,  and  7  is  the  ratio  of  specific  heats  of 


the  lasing  gas.  To  determine  the  relevant  characteristics  of  the  solutions  to  these  equations  one 
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Figure  4.1:  Geometry  of  laser 
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may  impose  boundary  conditions  p  =  0,  Q  =  0  at  the  walls  x  =  0,  x  =  d,  which  imply  that  p  and 
Q  may  be  expanded  in  discrete  Fourier  series  as 


p(x,t)  =  pn(t)  sin(fc„x) 

n 

Q(x,t)  =  ^  Qn(t)  sin (knx) 

n 


(4.2) 


where  kn  =  tx n/d  ( n  =  1,2,...)  are  the  allowed  spatial  frequencies.  These  are  not  the  most 
general  boundary  conditions,  but  the  pertinent  characteristics  of  the  solutions  will  not  depend  on 
the  boundary  conditions  except  as  discussed  specifically  below. 

Substituting  Eq.  (4.2)  into  Eq.  (4.1)  yields 


d_ 

dt 


\&_ 

dt 2 


+ ( kKu) 2 


Pn  =  -(7  -  l)k2nQn(t) 


(4.3) 


A  reasonable  approximation  is  to  consider  the  heating  rate  to  be  constant  over  the  interval 
0  <t  <  r,  where  r  is  the  pulsewidth  of  the  electrical  pumping.  Thus,  Q(x,t)  may  be  expressed 
as  H(x)/t ,  where  H{x)  is  the  total  electrical  single-pulse  energy  loading  per  unit  volume. 
Similarly,  Qn(t )  =  H„/t.  The  solution  of  Eq.  (4.3)  during  the  pulse  is  then 


pnit)  =  - 


uknt  —  sin(ufc„£) 
uknr 


(4.4) 


where  H  is  the  mean  energy  loading  per  unit  volume,  given  by 


H 


-If* 

dJo 


H(x)dx 


and 


hn  =  ~dJ0  (^)sin(fc*x)dr 


(4.5) 


(4.6) 


is  the  normalized  nth  Fourier  component  of  the  energy  loading  distribution. 

It  is  useful  to  consider  the  form  of  Eq.  (4.4)  in  the  short-time  and  long-time  regimes. 
kn  =  Trn/d  is  the  frequency  of  the  niU  acoustic  mode  of  the  channel.  For  each  mode  n  of 


wavelength  2n/kn  =  2 d/n,  one  may  make  the  following  two  observations:  (i)  When  not  enough 
time  has  passed  for  a  sound  wave  to  have  travelled  across  an  acoustic  wavelength  (uknt  <  1), 
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r^i 


t3.  This  is  the  well-known  transient  “t3"  regime  of  thermal 


sin  uknt  ~  uknt  so  that  pn 
blooming,  (ii)  For  longer  times  ( uknt  >>  1),  corresponding  to  many  acoustic  transits  across 
an  acoustic  wavelength,  the  blooming  reaches  the  steady  state,  and  pn  is  instead  proportional  to 
t.  Thus,  the  density  will  respond  more  strongly  to  small-scale  energy-loading  nonuniformities 
(with  uknr  >>  1)  than  to  larger-scale  nonuniformities  (with  uknr  <  1). 

In  the  following,  spatial  averages  (over  x)  will  be  denoted  by  <  . . .  >*,  and  time  averages 
by  <  ...  >t-  As  will  be  shown  in  the  next  section,  the  beam  quality  is  determined  by  the  spatial 
variance  of  the  density  perturbation,  averaged  over  the  pulse.  Using  Eq.  (4.4),  the  time-averaged 
variance  may  be  expanded  into  its  spatial  Fourier  components,  with  the  result 


«p‘>,>,=  \d‘  ]T  hli. 


where  B,  which  has  the  dimensions  of  a  density,  is  given  by 

(7-1)#' 


B  = 


U‘ 


(4.7) 


(4.8) 


and  In  =  1(9  =  uknr )  is  the  response  function  of  the  nth  acoustic  mode  of  the  channel  to  the 
energy-  loading  pulse  of  length  r: 


m  = 


6 3 


63  1  1 
—  -  2  sin  9  +  -9(  1  +  9  cos  9)  -  -  sin  29 


(4.9) 


The  behavior  of  I  as  a  function  of  the  normalized  pulsewidth  9  is  shown  in  Fig.  4.2.  It  is  seen 
that  the  transition  between  the  short-  and  long-time  regimes  is  relatively  sharp,  and  that  a  good 
approximation  for  the  function  1(9)  is 

Hr*  v.ti 

The  iA  dependence  at  short  times  corresponds  to  the  t&  dependence  of  p2  implicit  in  Eq.  (4.4), 
corrected  for  the  fact  that  at  a  constant  heating  rate  B2  would  be  proportional  to  t2. 

4.2  Beam  Quality 

The  energy  loading  nonuniformity  and  the  density  disturbance  which  it  causes  will  in  general 
not  be  independent  of  position  2  along  the  optic  axis  of  the  laser.  In  fact,  the  nonuniformities 
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NORMALIZED  PULSEWIDTH  0 


Figure  4.2:  Acoustic  response  versus  normalized  pulsewidth 
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will  occur  in  layers  along  z,  with  successive  layers  spatially  uncorrelated  with  each  other.  The 
phase  aberrations  from  each  layer  will  therefore  add  randomly  to  one  another,  much  like  a 
random  walk.  The  following  calculation  will  concentrate  first  on  the  effect  of  a  single  layer. 
The  addition  of  many  layers  will  be  carried  out  subsequently. 

The  phase  aberration  4>(x,t)  seen  by  a  ray  propagating  along  z  across  a  single  layer  is 
directly  proportional  to  p(x,t ): 

r  /3 

4>(x,  t)  =  ap(x,  t )  =  -71-  -  —  p(x,  t)  (4.11) 

A  Pa 

where  L0  is  the  correlation  length  along  z  (i.e.,  the  thickness  of  the  layer),  A  is  the  wavelength 
of  the  laser  light,  (3g  is  the  Gladstone-Dale  constant  (the  difference  between  the  refractive  index 
of  the  amplifying  medium  and  unity  at  standard  temperature  and  pressure  [STP]  ),  and  pA  is  the 
density  of  the  amplifying  medium  at  STP.  <j>(x,t)  may  be  expanded  in  a  Fourier  series  of  the 
form 

<t>(x,t)  =  '^2<f>n(t)sinknx  (4.12) 

n 

where  <f>n  =  opn,  and  the  spatial  variance  of  the  phase  aberration,  averaged  over  time,  may 
similarly  be  expanded  as 

<<  >x>t=  ~^>2J2hnIn  (4-13) 

n 

Here  $  =  aB  is  given  explicitly  by  the  following  convenient  expression: 


where  pA  is  the  standard  pressure  (1.01  x  105  Nt/m2,  or,  in  energy  density  units,  101  J/lit.)  Note 
that  $  depends  only  on  the  energy  loading  per  unit  volume  and  the  cell  length;  in  particular, 
it  is  independent  of  pressure.  However,  at  low  gas  pressures  and  correspondingly  low  energy 
loading,  $  can  be  reasonably  low.  A  typical  value  of  $  is  0.31  radians  at  L0  =  5  cm,  A  =  10 
^m,  /3g  =  3  x  10~4,  7  =  1.5,  and  H  =  10  J/lit. 

Beam  quality  is  determined  by  the  relative  value  of  the  peak  far-field  intensity,  known  as 
the  Strehl  ratio.5  In  principle,  the  Strehl  ratio  depends  on  both  the  phase  and  intensity  nonuni¬ 
formities  in  the  near  field.  In  practice,  however,  the  contribution  of  the  intensity  nonuniformity 
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to  deviations  of  the  Strehl  ratio  from  unity  is  usually  negligible.  For  the  case  at  hand  the  Strehl 
ratio  S(t)  may  then  be  written  as 


t 


S(t)  = 


<  >* 


(4.15) 


where,  as  above,  the  brackets  <  ...  >z  denote  a  spatial  average  over  the  near-field  beam  (i.e., 
over  x).  Expanding  the  exponent  in  powers  of  <j>  yields 


S{t)*l-(<<1>2  >x-  <<f>>l)  (4.16) 

to  order  <f>2.  In  general,  one  is  interested  in  the  temporal  average  of  the  beam  quality.  Let  S  be 
the  time  average  of  S(t).  It  follows  from  Eq.  (4.16)  that 


S  «  1-  <<  <f> 2  >i>«  +  «  <f>  >l>t 


(4.17) 


When  | <f>—  <  <f>  >  |  >  1  this  approximation  breaks  down,  and  higher  powers  of  <f>  must  be 
included.  If  the  energy-loading  distribution  satisfies  Gaussian-normal  statistics  then  the  power 
series  may  be  summed,  yielding 


S  «  exp  [-  <<  <f>2  >x>t  +  «<f>  >l>t 


(4.18) 


[For  other  statistical  distributions  of  the  energy  loading  Eq.  (4.18)  is  correct  only  to  order  <j>2; 
Eq.  (4.18)  and  all  subsequent  equations  describing  beam  quality  should  in  general  be  considered 
to  be  correct  only  to  this  order.] 

The  presence  of  hard  boundaries  confining  the  gas  in  a  fixed  volume  forces  the  mean  density, 
averaged  over  x,  to  be  a  constant,  independent  of  t.  Therefore  <  <f>  >z=  0.  Combining  Eqs. 
(4.10),  (4.13),  and  (4.14),  one  arrives  at  the  following  result: 


where 


S  «  exp 


/7TUT\4  V 

\~TJ  ^ 


n<  no 


3  d 


n0  = - 

7T  UT 
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is  approximately  the  number  of  acoustic  transits  across  the  channel  during  the  electrical  pulse. 
By  way  of  example,  for  d  =  2.4  cm,  r  =  30  //sec,  and  u  «  400  m/sec,  the  value  of  n0  is  ~  1.9. 

It  should  be  noted  that  beam  quality  is  often  measured  in  “number  of  times  diffraction 
limited”  ( XDL ).  The  relationship  between  XDL  and  S  is 

1 


XDL  = 


v/S 


(4.21) 


so  that  Eq.  (4.19)  implies  that 


XDL  «  exp  ^  +-$i 


\T)'  E  4  £  *1 } 

n<nQ  n>no  J  ) 


(4.22) 


The  above  results  apply  to  the  propagation  through  a  single  correlation  length  (along  z) 
of  the  density  disturbance.  Since  successive  layers  are  spatially  uncorrelated,  the  propagation 
through  several  correlation  lengths  results  in  the  addition  of  the  squares  of  the  phase  errors.  The 
generalizations  of  Eqs.  (4.19)  and  (4.22)  are  thus 


S  «  exp  { 

Z  jLq 


and 


XDL  «  exp 


’(=r)45>4*:+5£*il} 

n<«o  n>no  J  ) 

t  L  n<no  n>n0  J  ) 


(4.23) 


(4.24) 


where  L  is  the  total  distance  travelled  by  the  laser  beam  through  the  amplifier,  so  that  L/L0  is 
the  number  of  correlation  lengths  traversed  by  the  laser  beam. 

It  is  useful  to  break  the  right-hand  side  of  Eq.  (4.24)  into  two  parts,  as  follows: 


B 


XDL  «  exp  ^  +-$2— [6l  +  ds] 
L  o 


'} 


(4.25) 


where  hi  and  6s  are  the  contributions  from  large-  and  small-scale  nonuniformities,  respectively, 
in  the  electrical  loading: 

/itut\4 

n<r>o 

n>no 
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6l  is  the  contribution  of  transient  thermal  blooming  caused  by  large-scale  nonuniformities  (scale 
size  >  ut)  in  the  electrical  energy  loading.  Because  of  the  factor  of  (nurn/d)4  in  this  term,  it  is 
generally  small  (see  Fig.  4.2).  On  the  other  hand,  6s  is  due  to  the  steady-state  thermal  blooming 
caused  by  small-  scale  nonuniformities  (scale  size  <  ur)  in  the  energy  loading.  Because  of  the 
relatively  strong  steady-state  response  of  the  medium  to  the  smaller  scales,  it  is  65  which  must 
be  carefully  monitored  in  the  design  of  the  laser  amplifier. 

A  useful  alternative  expression  for  6s  is 

£>S  ~  ^  <  h2s  >x=  — [fts,  rmj]2  (4.27) 


where  /is(x)  is  the  small-scale  fractional  modulation  in  the  energy  loading: 


hs(x)  = 


H{x)  -  H 


H 


(4.29) 


Combining  Eq.  (4.28)  with  the  definition  of  $  yields 

XDL  «exp|+|  (y/*.1^—)2  ££o[A*™,l2j  (4.29) 

For  a  given  required  XDL,  Eq.  (4.29)  may  be  used  to  determine  the  acceptable  range  of  L0 
and  rms  variation  in  the  energy  loading.  By  way  of  example,  for  three  passes  through  a  125-cm 
amplifier  one  has  L  =  375  cm.  With  A  =  10 //m,  =  2  x  10-4,  7  =  1.5,  and  for  H  =  10  J/liter, 

Eq.  (4.29)  reduces  to 

XDL  *  exp  |+— o(.^Vm*)2  j  (4.30) 

where  L0  is  in  cm.  Eq.  (4.31)  implies  that  energy-  loading  nonuniformities  as  large  as  ~  20-30% 
rms  are  consistent  with  excellent  beam  quality  (1.1  XDL)  if  the  correlation  length  is  of  the  order 
of  10  cm  or  less. 

This  preliminary  analysis  thus  indicates  that  the  requirements  for  little  effect  on  beam  quality 
from  the  acoustics  driven  by  a  single  pulse  are  relatively  easy  to  achieve.  It  is  proposed  that 
a  more  detailed  analysis,  including  realistic  deviations  from  Gaussian-normal  statistics,  will  be 
carried  out  during  the  course  of  this  program. 
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4.3  Frequency  Noise 

A  second  important  effect  of  the  density  changes  during  the  pulse  is  changes  in  the  frequency 
of  the  amplified  laser  light.  In  general,  these  frequency  changes  may  be  divided  into  two  distinct 
parts:  chirp  and  frequency  “noise.”  The  chirp  uCh  is  the  change  in  the  instantaneous  central 
frequency  of  the  light,  while  the  noise  8v  is  associated  with  deviations  from  the  instantaneous 
central  frequency,  i.e.,  creation  of  Fourier  components  at  other  frequencies.  For  a  laser  amplifier 
with  solid  walls,  the  chirp  should  be  nearly  zero  if  the  laser  beam  uniformly  fills  the  entire  region 
between  the  walls  (i.e.,  if  the  laser  intensity  is  independent  of  x).  Deviations  from  uniformity 
may  produce  some  chirp. 

4.4  Mode-Medium  Instability 

Another  effect  that  could  in  principle  prove  relevant  to  controlling  the  beam  quality  of  a 
pulsed  C02  laser  is  the  mode-medium  instability  (MMI),  which  has  been  described  in  detail  in 
Ref.  4.  The  MMI  is  due  to  the  interaction  between  the  laser  intensity  distribution,  the  dependence 
of  the  heating  rate  on  the  laser  inversion,  the  acousto-density  response  of  the  gaseous  medium 
to  the  heating  distribution,  and  the  refraction/diffraction  of  light  by  the  density  disturbance. 

It  is  anticipated  that  a  C02  laser  amplifier  of  the  type  used  for  laser  radar  will  not  suffer 
any  significant  loss  from  MMI.  As  was  shown  by  Flusberg  et.  al.,  the  instability  is  important  in 
a  laser  oscillator,4  but  does  not  cause  a  loss  of  beam  quality  in  a  laser  amplifier.6,7  In  addition, 
even  in  an  oscillator,  the  rate  at  which  the  instability  grows  depends  on  density;  at  a  density 
of  30-40  torr  the  degradation  in  oscillator  beam  quality  should  occur  on  a  time  scale  of  ~  100 
^sec,  hence  there  should  be  no  significant  effect  on  the  beam  quality  of  a  30-/*sec  pulse. 

4.5  Persistent  Acoustic  Disturbances 

Another  important  effect  on  beam  quality  is  that  of  persistent  acoustic  disturbances,  which 
are  created  by  the  heat  deposition  during  each  pulse.  Between  pulses  these  disturbances  relax 
from  the  combined  effects  of  viscous  damping  in  the  gas  and  imperfect  acoustic  reflections  from 
the  solid  boundaries.  For  a  repetitively  pulsed  laser,  each  successive  pulse  “feeds”  the  various 
acoustic  modes  of  the  laser  cavity.  Each  mode  decays  between  pulses  but  is  driven  again  by 
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successive  pulses.  The  result  is  a  steady  state  between  the  driving  impulses  due  to  the  pulsed 
heat  deposition  on  the  one  hand  and  the  damping  by  viscosity  and  finite  reflectivity  on  the  other. 
For  a  small  enough  wall  reflectivity  the  acoustics  should  damp  away  so  rapidly  that  no  significant 
buildup  should  take  place. 
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